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ABSTRACT 

In the context of star formation through fragmentation of an extremely metal-deficient protogalactic 
cloud, the gravitational collapse of filamentary gas clouds is explored with one-dimensional numerical 
hydrodynamics coupled with non-equilibrium chemistry of H2 and HD. It is found that the cloud evolution 
is governed mainly by the initial central density (n^o) and H2 abundance (a:H2,o)- In particular, the 
evolution of low-density filaments (tt-co ^ 10^ cm~'^) bifurcates at a threshold H2 abundance of 2;H2,cr — 
3 X 10^^, beyond which HD cooling overwhelms H2 cooling. The contraction of a filament with nc,o ^ 10^ 
cm~'^ and Xh2,o ^ XH2,cr is strongly decelerated when the central density (ric) reaches a critical density 
of HD at which LTE level populations are achieved, and therefore the filament is expected to fragment 
at ~ 10^ cm^"^. The fragment mass is assessed to be « IOMq. In contrast, the contraction of a filament 
with ricfl < 10^ cm~^ and a;H2,o ^ 2;H2,cr is regulated by H2 cooling. In this case, the filament tends 
to fragment at lower density as ~ 10'' cm~'^ owing to the low critical density of H2, and the fragment 
mass is as high as « 10^ Mq. For a high-density filament with n^o ^ 10^ cm~'^, the temperature stays 
at a relatively high value because both H2 and HD cooling saturate, and the cloud evolution is governed 
by H2 cooling. The contraction of a high-density filament is accelerated by effective three-body H2 
formation when the density reaches 10^~^ cm~'^. The fragmentation is not expected to take place until 
the cloud becomes opaque in H2 lines at n^o ~ 10^^~^^ cm~^, so that the fragment mass is reduced to 
1 — 2 Mq. As a result, the stellar initial mass function (IMF) could be bimodal and deficient in sub-solar 
mass stars, where the high mass peak is around lOM© or IO^Mq, dependently on n^o and xb_2,o- If 
the protogalactic clouds are ionized by UV radiation or strong shocks, the H2 abundance could exceed 
xn^ cr — 3 X 10~^ by reactions of H -I- e — * H^ + hi' and H + H^ ^ H2 + e. The high mass peak would 
then be O(10)Mo. 

Subject headings: cosmology: theory — galaxies: formation — hydrodynamics — ISM: clouds — stars: 
formation 



1. INTRODUCTION 

The first generation of stars, Population HI stars, are 
thought to form from almost metal-free gas with metallic- 
ity of ~ 10-^°Zq (see e.g., Carr, Bond, & Arnett 1984; 
Carr 1994). Such Population III stars may have produced 
heavy elements at redshifts z > 10. Recent observations of 
quasar absorption spectra show that the metal abundance 
in intergalactic space is « 10~^ Zq (e.g., Cowie & Songaila 
1998). Also, it is widely accepted that a significant por- 
tion of old halo stars in our Galaxy have metallicity lower 
than 10~^Zq (Beers, Preston, & Shectman 1992). In ad- 
dition, some blue compact dwarf galaxies are known to be 
extremely metal-poor (« 10~^Zq, Kunth & Sargent 1986; 
Pustilnik et al. 2001). Thus, star formation in protogalax- 
ies must have proceeded in very metal-deficient environ- 
ments. 

When the metallicity is lower than Z < 10^^ Zq, the 
cooling by heavy elements is less effective than cooling 
by hydrogen and helium (e.g., Yoshii & Sabano 1980; 
Bohringer & Hensler 1989; Omukai 2000; Nishi & Tashiro 
2000). For the formation of Population HI stars from 
such metal-deficient gas, the cooling by primordial hy- 
drogen molecules (H2) plays an essential role (Matsuda, 



Sato, & Takeda 1969; Yoneyama 1972; Hutchins 1976; 
Silk 1977; Yoshu & Sabano 1980; Carlberg 1981; Lepp 
& ShuU 1984; Palla, Salpeter, & Stabler 1983; Yoshii & 
Saio 1986; Shapiro & Kang 1987; Uehara et al. 1996; 
Haiman, Thoul & Loeb 1996; Nishi et al. 1998; Abel et 
al. 1998; Bromm, Coppi, & Larson 1999; Abel, Bryan, & 
Norman 2000; Nakamura & Umemura 1999a, 2001, here- 
after Papers I and II). Recent studies have revealed that 
star formation in primordial gas is considerably different 
from present-day star formation, implying that the stellar 
initial mass function (IMF) might be different from the 
Salpeter-like IMF and the IMF might be time-varying in 
the course of galaxy evolution from the early collapsing 
stages to the present day (see also Larson 1998; Zepf & 
Silk 1996; Chabrier 1999). 

Besides hydrogen molecules, deuterated hydrogen 
molecules (HD) can be a significant coolant (e.g., Galli 
& Palla 1998). Although HD is less abundant than H2 
([HD/H2] ~ 10^3 _ 10-4)^ HD has a finite dipole mo- 
ment and thus higher radiative transition probabilities 
than H2. For example, the lowest rotational transitions 
of H2 and HD have radiative transition probabilities of 
A20 = 3 X 10^^^ s~^ and Aio ~ 5x 10^^ s^^, respectively. 
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Also, the corresponding excitation energies of H2 and HD 
are AE2o/k = 510 K and AEio/k = 128 K, respectively, 
where k is the Boltzmann coefRcient. Thus, HD can lower 
the gas temperature down to T < lOOK (e.g., Puy & Sig- 
norc 1996; Bouglcux & GalU 1997; GalU & Palla 1998; 
Flower et al. 2000). 

The first pregalactic objects should have collapsed at 
redshifts of ~ 10 - 10^ and have masses of 10^ - 10^ Mq 
in a cold dark matter cosmology (Tegmark et al. 1998; 
Fuller & Couchman 2000). When such first objects arc 
virialized and the gas temperature ascends to 10^ — lO"' K, 
H2 formation is promoted and the fractional abundance is 
raised to = 10~^ — 10^"^. However, for this H2 abun- 
dance, the H2 cooling cannot lower the temperature down 
to 100 K. Thus, the cloud evolution is basically controlled 
by H2 cooling and HD cooling is not important (Nakamura 
& Umemura 1999b, 2000; see also Lepp & Shull 1984). 

But, in some situations, the H2 abundance can be raised 
further. Ferrara (1998) argued that in dense shells formed 
behind supernova shocks, H2 molecules form efficiently ow- 
ing to reduced recombination of free electrons (see also 
Shapiro & Kang 1987). In the dense shells, the H2 concen- 
tration reaches 6 x 10"'^, and consequently the gas temper- 
ature descends to ~ 100 K. At such low temperature, HD is 
a more efficient coolant than H2. The primordial star for- 
mation regulated by HD cooling in shocked shells is stud- 
ied by several authors (Uehara & Inutsuka 2000; Machida, 
Fujimoto, & Nakamura 2001). Another possibility is star 
formation from photoionized gas. Corbelli, Galli, & Palla 
(1998) studied the effects of UV background radiation on 
the thermal evolution of the protogalaxics (sec also Mur- 
ray & Lin 1989; Haiman, Rees, & Loeb 1996). They found 
that at redshifts lower than z ~ 1 — 2, protogalactic clouds 
are self-shielded against the UV background radiation, and 
H2 formation is promoted up to Xu^ ~ 10~^ with the help 
of abundant free electrons due to the UV radiation. Even 
at higher redshifts {z > 2), if pregalactic clouds more mas- 
sive than 10^^ M© collapse under a UV background, they 
can be eventually self-shielded and abundant H2 (> 10~^) 
forms, so that the temperature descends down to ~ lOOK 
(Susa & Umemura 2000). Hence, HD cooling is likely to 
be important for the star formation in the protogalactic 
systems. In this paper, we extensively examine the effects 
of the HD cooling on the formation of stars in protogalax- 
ies, and elucidate the role of HD molecules for the stellar 
initial mass function (IMF) there. 

This paper is organized as follows. In §2, we describe 
model and numerical methods which are basically the same 
as those of Paper II. Numerical results are presented in 
§3. In §4, implications for the IMF of Population III and 
metal-deficient stars are discussed. 

2. MODEL AND NUMERICAL METHODS 

In Papers I and II, we studied the fragmentation of fil- 
amentary clouds by incorporating the H2 cooling. In this 
paper, we examine the effects of HD cooling on the evo- 
lution of the filaments. Our numerical model and method 
are the same as those of Paper II except for the inclusion 
of deuterium chemistry. In the following, we briefly review 
our numerical model. See §2 of Paper II for more detail. 

We deal with the following 14 species: e, H, H+, H~, H2, 
H^, He, He+, He++, D, D+, D", HD, and HD+. The mass 



fraction of He is set to 0.24 of the total mass. The D abun- 
dance is set to 4 X 10~^ by number, which is consistent with 
recent observations of the deuterium Ly-a feature in the 
absorption spectra of high-redshift quasars (e.g., Tytler et 
al. 1996; O'Mcara et al. 2000). The reaction rate coef- 
ficients for the deuterium chemistry are given in Table 1, 
while the rate coefficients for other species are the same as 
those of Paper II. 

We take into account the following thermal processes: 
(1) H cooling by radiative recombination, coUisional ion- 
ization, and collisional excitation, (2) H2 line cooling by 
rotational and vibrational transitions, (3) cooling by H2 
collisional dissociation, (4) heating by H2 formation, and 
(5) HD line cooling by rotational transitions. 

The HD line cooling is computed by the same method as 
the H2 line cooling described in Paper II, which includes 
the escape probability method in an optically-thick regime. 
For the collisional dcexcitation rates for HD. we consider 
both H-HD and H2-HD collisions, using analytical fits of 
Galh & Palla (1998) and Flower & Roueff (1999). The 
first 6 rotational levels are taken into account. Recently, 
Flower et al. (2000) derived an updated HD cooling func- 
tion in an optical thin regime, which includes up to 8 rota- 
tional levels. We have checked the difference by employing 
their cooling function for several models described in §3, 
and found that the resultant temperatures agree with the 
present results to within ~4%. 

We consider an infinitely long cylindrical gas cloud 
which collapses in the radial direction. The initial tem- 
perature is assumed to be spatially constant. The initial 
abundances of electron and D"*" are set to a:;e = 5 x 10~^ 
and xy)+ = 1 x 10^^, respectively. (The numerical results 
do not depend sensitively upon .t,. and . See also §2 
of Papers I and II.) At the initial state, the relative abun- 
dances of H-, H+, He+, He++, D", HD, and HD+ are set 
to zero for simplicity. The abundances of other species are 
determined by the conservation of mass and charge. The 
initial density profile in the radial direction is assumed to 
be 

Po = Pc,o{l + ryRiy\ (1) 

where Rq = -y/2/A;To/(7rGpc,oM) is the effective radius, Pc,o 
is the central mass density, Tq is the initial gas tempera- 
ture, fi is the mean molecular weight, and / is the ratio of 
the gravitational force to the pressure force. When / = 1, 
the density distribution agrees with that of an isothermal 
cylinder in hydrostatic equilibrium. When the filamentary 
cloud forms through gravitational fragmentation of a par- 
ent sheet-like cloud, / is expected to be w 2 (see §2.2 of 
Paper I). The radial infall velocity is given by 

Vrfi = -vor (^Ro + \[^+r^^ , (2) 

where vq is constant and is set to the initial sound speed. 
We neglect the effects of dark matter because after virial- 
ization of the parent system, the local density of baryonic 
gas is likely to be higher than the background dark matter 
density owing to radiative cooling (e.g., Umemura 1993). 
The model is thus specified by four parameters: Uc^, Tq, 
/, and the initial H2 fractional abundance x^^fi- As stated 
in §1, high H2 abundance as xh2,o ^ 10~^ can be achieved 
if a parent cloud undergoes ionization prior to the collapse. 
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In Paper II, we followed the collapse and fragmentation 
of the primordial filaments with one-dimensional and two- 
dimensional axisymmetric simulations. The numerical re- 
sults of Paper II showed that the fragment masses derived 
from the two-dimensional simulations are in good agree- 
ment with the estimate based on one-dimensional simula- 
tions. Therefore, in this paper, we pursue one- dimensional 
simulations on the collapse of the filaments to estimate the 
fragment masses. 

We calculated more than 3000 models by choos- 
ing the model parameters nc,o, ^o, /, and xh2,o as 
logio(nc,o/cm-3) = 1.00, 1.33,1.67, 2.00, 6.00, Tq = 
200, 300, 400 K, f = 1.5, 2.0, 2.5, 6.0, and xr^ o = 
1 X 10-^ 3 X 10-^ 6 X 10-4, 1 X 10-3, . . . , 1 X 10-2, respec- 
tively. 

3. NUMERICAL RESULTS 

In Papers I and II, we studied the fragmentation of the 
primordial filaments incorporating H2 cooling, and it was 
shown that the fragmentation is bifurcated by a threshold 

initial density of « lO^cm^'^, below which as massive stars 
as « IO'^Mq form, and above which less massive stars with 
w IM0 form. 

In this section, we reexamine the collapse of the fila- 
ments, including the HD cooling as well as H2 cooling. 
As shown below, there is a threshold initial H2 concentra- 
tion, above which HD cooling predominantly regulates the 
cloud evolution. We also show that the HD cooling does 
not play an important role for high-density gas. Thus, the 
evolution of the filaments is classified into three cases, de- 
pending upon the initial density, ^co, and initial H2 abun- 
dance, Xu2fi', (1) low-density filaments with high XHa.o, (2) 
low-density filaments with low a;H2,05 and (3) high-density 
filaments. 

In the following, we first compare the HD cooling and 
H2 cooling on the ricfi — xn^fi diagram to clarify the impor- 
tance of HD cooling. Next, we show the numerical results 
of one-dimensional simulations. In the models shown in 
this section, the initial temperatures are set to be equi- 
librium temperatures of the filaments which decrease with 
increasing a;H2,o- 

3.1. Threshold H2 Abundance 

HD cooling depends sensitively on the initial cloud den- 
sity and H2 abundance. To clarify the effect of HD cool- 
ing, we compare the HD cooling rate (Ahd) and the H2 
cooling rate (Ahj) in Figure 1. The abscissa and ordinate 
indicate the initial gas density and the H2 abundance, re- 
spectively. The solid lines show the contour curves of the 
HD-to-H2 cooling ratio (Ahd /Ana)- When evaluating the 
cooling rates, the gas temperatures are determined itera- 
tively to satisfy the condition of tcoo\ = ifrag; where tcooi 
is cooling time and tfrag is fragmentation time as defined 
in §3.2. The resultant temperatures are also shown by 
dashed lines. The HD abundance is taken to be propor- 
tional to the H2 abundance as Xud/xh2 = 1 x 10-^, which 
is consistent with those in the models shown in §3.3 and 
3.4. This figure shows that the contribution of HD cooling 
rate is largest around ric.o — 10^ cm^"^, which is almost 
comparable to the critical density of HD, beyond which 
the rotational level populations achieve the LTE. 

If the cloud has the initial H2 abundance lower than a 



threshold value of XHa.cr = 3 x 10""^, HD cooling does not 
play a significant role in the thermal evolution of the cloud 
because Ahd < Ahj • The cloud temperature is then equal 
to the equilibrium temperature shown by the dashed lines. 
As shown later, the H2 abundance in a collapsing filament 
stays almost constant before the three-body reactions pro- 
ceed at 10®-^ cm-3. Thus, if an evolutionary path of such 
a cloud is superimposed on Figure 1, then it is almost par- 
allel to the abscissa. For the filaments with n < 10^ cm-^ 
and .THa ^ •^H2,cr, the evolutionary path enters into the re- 
gion of Ahd > Ah2 in the course of contraction even if Ah2 
is larger than Ahd at the initial state. In actual evolution, 
as shown below, the pressure force retards the contraction 
of the filament. In particular, when the density exceeds 
HD critical density, nnD.cr ^ 10*^''' cm^'', the contraction 
slows and the actual cloud temperature stays below 100 
K. Consequently, HD cooling continues to control the con- 
traction until the cloud becomes opaque to the HD lines. 
On the other hand, for the filaments with nc.o ^ cm^"^ 
and a;H2,o ^ 2^H2,cr; the evolutionary path goes into the re- 
gion of Ahd < Ah2 even if Ahd is larger than Ah2 at the 
initial state. As a result, for the initial condition in the 
gray region in Figure 1 , the HD cooling plays an important 
role in the thermal evolution of the gas. 

3.2. Low-Density Filaments with high H2 abundance 
(nc,Q < lO-'^cm-^ and .th2,o'> 3 x lO"-"*; 

As a typical example in which the HD cooling controls 
the contraction, we show the evolution of the model with 
nc,o = 10 cm-3. To = 200 K, / = 2, and a;H2,o = 3 x 10-^. 

Figures 2a, 2b, 2c, and 2d show the time evolution of (a) 
the temperature, (b) H2 {solid line) and HD (dashed line) 
abundances, (c) cooling rates by H2 {dashed line) and by 
HD {dotted line), and the total {solid line), (d) the cooling 
timescale {solid line), the timescale of dynamical contrac- 
tion {dashed line), and the fragmentation timescale {dot- 
ted line), respectively, as a function of the central density 
[Uc = pc/{jimH)]- Since the central density monotonically 
increases with time, the abscissa corresponds to the evo- 
lution time. The contraction and cooling timescale are, 
respectively, defined as idyn = p/p and tcooi = Et/{At), 
where Et and At are total internal energy and net cooling 
rate, respectively. The fragmentation time, ffrag, is defined 
as 2.5 times the inverse of the growth rate of the fastest- 
growing mode (eq [38] of Nakamura, Hanawa, & Nakano 
1993), i.e., ifrag = 5.17(27rGp)"i/^ Note that this defini- 
tion of tfrag is 2.5 times larger than that of Paper II. In 
practice, a filament does not fragment in one- dimensional 
calculations. Therefore, tfrag should be regarded as a mea- 
sure of the fragmentation epoch as well as the free-fall 
time-scale. 

At the early stages of the contraction, H2 cooling dom- 
inates HD cooling, and the equilibrium temperature is 
lowered to T ~ 100 K owing to the high H2 abun- 
dance. When the density reaches the critical density of H2, 
?iH2,cr ~ 10"^-** cm-"^, the H2 cooling is saturated and then 
HD cooling becomes dominant as the density increases. 
When the HD cooling becomes effective, the equilibrium 
temperature descends to ^50 K owing to the lower exci- 
tation temperatures of HD. Thereafter, the temperature 
stays nearly constant at ~50 K during the contraction. 
When the density reaches 

^HD.cr 10 cm , the cloud 
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contraction tends to become quasistatic by the combina- 
tion effect of (1) tlie dynamical stability of a cylindrical 
cloud and (2) the characteristic of line cooling. First, a 
polytropic cylinder in hydrostatic equilibrium with P (x 
is gravitationally stable in the radial directions, if the ef- 
fective 7 is greater than unity. This means that an isother- 
mal cylinder is marginally stable. Second, at lower den- 
sity than the critical density, the line cooling rate is nearly 
proportional to because the collisional excitation rate 
balances with spontaneous emission rate, while at higher 
density than the critical density, the cooling rate is pro- 
portional to n because the level populations reach LTE. 
Hence, tcooi is inversely proportional to the density for 
n < "-HD.cr, while tcooi is constant for n > riHD.cr if the 
temperature is constant. By these effects, the compres- 
sional heating strongly brakes the dynamic contraction 
when n > riHD.cr and then the cloud evolves quasistati- 
cally. 

When the density reaches ~ 10^ cm^'^, the cloud be- 
comes optically thick in the HD lines and the quasistatic 
contraction almost ceases. However, before the central 
density reaches this density, the filament is likely to frag- 
ment. As shown by two-dimensional calculations in Pa- 
per II, it is anticipated that the fragmentation takes place 
when the dynamical timescale becomes a few times longer 
than the fragmentation timescale (see Paper II for more 
detail). From Figure 2d, this cloud is expected to undergo 
the fragmentation when the density reaches 10^ cm^''. 

Uehara & Inutsuka (2000) estimated the minimum frag- 
ment mass when the HD lines become opaque. However, in 
their model, the cloud becomes optically thick to the HD 
lines when the density reaches « 3 x 10^° cm^'^, which is 
about 30 times higher than the present result. To under- 
stand this discrepancy, we compare our evaluation of the 
optical depth with the LVG approximation (Goldreich & 
Kwan 1974) with some modification. The optical depth of 
the LVG approximation is originally given by 

hcBj^j+mj f gjnj+i\ 

where h is the Plank constant, c is the speed of light, 
Bj j^i is Einstein's B coefficient of the transition J 
J + 1, \dv/ds\ is the velocity gradient, nj and gj are the 
number density and the Gaunt factor of the J level, respec- 
tively. However, in the present calculations, the velocity 
gradient is small compared to the thermal width and the 
gradient of the thermal velocity is dominant. Hence, we re- 
place the velocity gradient by avth/Rj = ay^nGpc, where 
a, Uth, and Rj are a nondimensional numerical constant, 
the thermal velocity, and the filament radius, respectively. 
If we set the temperature to 60 K, the optical depth of 
the most effective cooling line ( J = 1 — > 0) becomes unity 
when the density reaches 1 — 3 x 10^ cm"'^ for a = 1 — 2. 
This is consistent with the present result. Also, it should 
be noted that our definition of the fragment mass is dif- 
ferent from that of Uehara & Inutsuka (2000). We adopt 
the mass contained within one wavelength of the fastest- 
growing linear perturbation which is consistent with our 
two-dimensional numerical results (sec §4), while Uehara 
& Inutsuka (2000) adopted the Jeans mass which is about 
ten times smaller than our definition at the same density 
and temperature. As a result, their estimation of the frag- 
ment mass is by about two orders of magnitude smaller 



than our estimate if we assume the cloud to fragment at 
the stage optically thick to the HD lines. However, in prac- 
tice, the cloud is likely to fragment at the earlier stages 
around the critical density of HD as shown above. 

The evolution of the other models with different initial 
parameters (10 cm'^ < ric o < 10"* cm'^, 200 K < Tq < 
400 K, 1.5 < / < 6, and 3 x IQ-^ < xu.fi < lO"") is 
qualitatively similar to that of the model shown in Figure 
2. When the initial H2 abundance is as high as 10~^, the 
HD cooling is more effective; than H2 cooling at the initial 
state. Thus, the temperature goes down to ~ 50 K at the 
very early stages of the evolution. 

3.3. Low-Density Filaments with Low H2 Abundances 
Ca:H.,o<3x IQ-^) 

As a typical example of this case, we show the evolution 
of the model with nc,o = 10 cm'^, Tq = 400 K, / = 2, and 
2;h2,o = 1 X 10^* in Figure 3. The notation of the figure 
is the same as that of Figure 2. 

In this model, Ahd < Ah2 during the contraction, and 
thus HD cooling does not play a significant role in the ther- 
mal and dynamical evolution. The temperature remains 
between 100 and 500 K over 11 orders of magnitude in den- 
sity until the H2 lines become optically thick. During the 
contraction, the HD abundance is nearly proportional to 
that of H2 (xhd ^ lO^^'aJHa) because HD molecules form 
primarily via the reaction D + H2 — > H + HD. Thus, when 
the three-body H2 formation becomes effective, almost all 
the D atoms are processed into HD. It should be noted 
that the above reaction is sensitive to the tcmpcratTirc, and 
for lower temperature, the reactions D+ + H2 ^ H+ + HD 
and HD+ -|- H ^ D+ + HD are more effective. 

The contraction proceeds quasistatically after the den- 
sity reaches the critical density of H2 (riHs.cr = 10^~^ 
cm~^). When the three-body H2 formation becomes ef- 
fective, the contraction is accelerated again and then de- 
celerated when the cloud becomes optically thick to the H2 
lines (ric ~ 10^^ cm~^). The cloud is expected to fragment 
after the density reaches the critical density of H2. 

The evolution of the other models with different param- 
eters (10 cm-3 < nc,o < 10^ cm'^, 200 K < T < 400 K, 
1-5 < / < 6, and lO"-* < xh^.o < 10"^) is qualitatively 
similar to that of this model. 

3.4. High-Density Filaments 

As a typical example of this case, we show the evolution 
of the model with nc,o = 10^ cm'^, Tq = 200 K, / = 6, 
and xii2,o = 1 x 10~^ in Figure 4. The notation of the 
figure is the same as that of Figure 2. 

Although the initial H2 abundance is as high as the 
threshold H2 abundance 3 x 10^'^, the evolution is essen- 
tially the same as those of the models with low xh^^q (see 
§3.2 of Paper II). S ince tIq q is already higher than tih^ qj- 
and nHD.cr, the temperature increases as the collapse pro- 
ceeds. Therefore, even if HD cooling is more effective than 
H2 cooling at the initial state, H2 cooling becomes domi- 
nant as the collapse proceeds. In other words, the evolu- 
tion is essentially the same as that without HD cooling. 

The contraction time does not become longer than the 
fragmentation time until the H2 lines become optically 
thick (ric ^ 10^^ cm~^) due to the rapid increase in the H2 
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abundance via the effective three-body reactions. There- 
fore, the fragmentation is not expected to take place until 
the cloud becomes opaque to the H2 lines. 

The evolution of the other models with different param- 
eters (ric > 10^ cm-^, 200 K < T < 400 K, 3 < / < 6, 
and 10-'' < 2^H2,o ^ 10 ^) is qualitatively similar to that 
of this model. 

It is noted that for the high-density filaments with low 
f (e.g., / < 2 for To = 300 K), t^yn becomes comparable 
to or longer than th^g before ric roaches 10** cm^'^. This is 
because when the parameter / is small (or the initial line 
mass is small), the equilibrium line mass of the filament, 
which is proportional to the temperature, approaches the 
total line mass of the cloud before the three-body reactions 
become significant. The contraction of such a cloud is de- 
celerated significantly around Uc ~ 10* cm"^, and thus 
the cloud is likely to fragment around this density. 

4. DEPENDENCE OF THE FRAGMENT MASS ON THE 
INITIAL MODEL PARAMETERS 

In this section, we estimate the fragment masses with 

our numerical results. As shown above, there is a threshold 
H2 abundance, beyond which the HD cooling can play a 
key role in the thermal evolution of primordial gas clouds, 
and also there is a threshold initial density, below which 
the filament would fragment at the critical density of H2 
or HD. Hence, as anticipated from the above numerical 
results, the fragment mass is also divided into three cases, 
depending upon the initial H2 abundance and density. 

Comparing one-dimensional results with the two- 
dimensional results shown in Paper II, we found that the 
fragmentation takes place when t^yn reaches 2—3 times the 
inverse of the growth rate of the fastest-growing perturba- 
tion. We thus assume that a filament fragments into dense 
cores when the dynamical time reaches the fragmentation 
time, idyn = ifrag- The fragment masses are computed 
as A'/frag ^ rricffX, where meg = J 2'Krpdr is a line mass 
which is obtained by integrating the density from the cen- 
ter to the radius at which the density takes O.lric and A is 
a longitudinal wavelength of a fastest-growing linear per- 
turbation (cq.[38] of Nakamura et al. 1993). This mass is 
about ten times larger than the Jeans mass. 

Figures 5a and 5b show the fragment mass for the mod- 
els with xhs.o = 1 X 10"^ and 3 x 10"^, with assuming 
equilibrium initial temperature. The abscissa and ordi- 
nate denote the initial central density and the parame- 
ter /, respectively. The solid lines denote the contours of 
the fragment mass with adjacent numbers in units of Mq. 
The thick dashed lines show the boundary at which the 
HD cooling rate is equal to H2 cooling rate at the epoch 
of fragmentation. In the regions to the left of the dashed 
lines, HD cooling is more effective than H2 cooling. Fila- 
mentary clouds are expected to form by the fragmentation 
of sclf-gravitating disks. According to numerical simula- 
tions, when the filaments form by the growth of the most 
unstable modes, their line masses tend to be twice the line 
mass for hydrostatic equilibrium, /Zeq = 2Cg/G, where Cg 
is a sound speed. For comparison, the line masses of the 
filaments with 1.5/icq, 2/icq, and 2.5/icq are shown by dot- 
dashed lines. The line mass of the equilibrium filament 
depends only on the equilibrium temperature that is de- 
termined by the balance between the cooling time and the 



fragmentation time of the self-gravitating disk. Here, the 
fragmentation time is defined as the inverse of the most 
unstable linear perturbation (see Larson 1985). The to- 
tal line masses of the filaments are proportional to /Tq. 
Therefore, the ordinate in Figure 5 corresponds to the to- 
tal line masses of the filaments. 

For the models with xhs.o < a;H2,cr ~ 3 x 10~^ (Fig. 
5a), the distribution of the fragment mass is quite similar 
to the case without HD (Paper II), because HD cooling 
docs not play an important role in the thermal evolution 
of the filaments. For higher initial density and/or larger 
/, the fragment mass is lower. The maximum and min- 
imum masses are estimated as ^ IO^Mq and 1 — 2Mq, 
respectively. The former corresponds to the Jeans mass 
at the stage at which the density reaches a critical density 
of H2. The latter corresponds to the Jeans mass at the 
stage at which the cloud becomes opaque to the H2 lines. 
[The maximum mass is a few times smaller than that of 
Paper II because of the higher initial H2 abundance. For 
the models with Xhs.o = 1 x 10^"'', the mass distribution 
as well as the maximum mass is in good agreement with 
that of Paper II. See Figure la of Nakamura & Umemura 
(2002).] 

There is a steep boundary at nc,o ~ 10** — 10^ cm~"^ 
in the distribution of the fragment mass for / > 3. For 
the models with nc,o ^ 10^ cm~^, the fragment masses 
take their minima at 1 ~ 2 M©, whereas, for the mod- 
els with nc.o ^ 10'^ cm^'^, they are greater than ^ 10^ 
M0 . This sensitivity in the fragment mass comes from the 
rapid increase in the H2 abundance due to the three-body 
reactions. For the models with low densities (nc,o ^ 10^ 
cm^"^), the contraction becomes quasistatic when the den- 
sity reaches the critical density of H2 and then linear den- 
sity fluctuations can grow nonlinearly before the three- 
body H2 formation becomes dominant (nc,o ^ 10*~^ 
cm^'"*). In contrast, for the models with high densities 
(nc,o ^ 10^~^ cm"*^ and / ^ 3), the contraction time does 
not exceed the fragmentation time until the H2 lines be- 
come optically thick at n 10^^ — 10*^ cm~^. 

On the other hand, when the initial H2 abundance is 
higher than 3 x 10~^ (Fig. 5b), the HD cooling is more 
effective than the H2 cooling for low-density filaments. 
Thus, the maximum mass of the low-density region is 
reduced to a few tens Mq (The fragment mass depends 
weakly upon the initial H2 abundance. For example, for 
the models with xhs.o = 1 x 10^^, the fragment mass is 
around 10 Mq in the low-density region.). The reduced 
maximum mass is then related to the Jeans mass at the 
stage at which the density reaches a critical density of 
HD. The minimum mass does not change because HD is 
not a dominant coolant in the evolution of dense filaments. 
Therefore, similarly to the models with low H2 abundance, 
the dependence of the fragment mass on the initial density 
exhibits a steep boundary around nc,o = 10* — 10^ cm~^ 
for / > 3. 

As mentioned in §3.4, for the high density filaments with 
low /, the contraction time temporally becomes longer 
than the fragmentation time just before the three-body 
reactions become important. Thus, fragmentation is ex- 
pected to take place before the three-body H2 formation 
becomes efficient. As a result, the fragment masses are 
relatively large (the lower-right regions in Figure 5). 
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In summary, the distributions of the fragment mass has 
a steep boundary around nc,o ^ 10^~^ cm~^ for / > 3, ir- 
respectively of the initial H2 abundance. This implies that 
the IMFs in very metal-deficient gas are likely to be bi- 
modal if both low-density and high-density filaments form 
stars. The low-mass peak is around a few Mq, which is not 
sensitive to the abundance of H2 formed in a parent cloud. 
The high-mass peak is « lO^M© if ajRa.o ^ 3 x 10~^, while 
it is « IOM0 if XH2,o > 3 X 10-3. 

5. DISCUSSION 

First, we briefly discuss the role of HD for the formation 
of the very first stars (Pop III stars). In the bottom- up 
scenarios like cold dark matter models, when the first pre- 
galactic objects with masses of 10^ — 10^ M© are virialized 
at redshifts of 2; ~ 10 — 10^, the H2 abundance reaches 
at most 10—* — 10"^ which is lower than the threshold 
II2 abundance. Therefore, the cloud evolution is basi- 
cally determined by H2 cooling rather than HD cooling 
(Nakamura & Umemura 1999b, 2000). If the primordial 
D abundance is a few times as high as the value we as- 
sumed (a;D = 4 X 10^'"'), HD cooling may play a role in 
the thermal evolution of the pregalactic clouds. However, 
recent observations of quasar absorption spectra find an 
observed primordial D abundance as low as 3 — 4 x 10"^ 
(Tytler ct al. 1996), implying that HD cooling is not im- 
portant for the first star formation. 

Next, wc give some implications for star formation in 
metal-deficient primordial galaxies. Recently, Susa & 
Umemura (2000) investigated the pancake collapse of pre- 
galactic clouds under UV background radiation. They 
found that once the pancake disk is shielded against exter- 
nal UV radiation in the course of contraction, H2 molecules 
form efficiently via the H~ reaction with abundant free 
electrons produced by UV background, and the resultant 
abundance reaches .TH2 ~ 3 x 10^'' (see also Shapiro & 
Kang 1986). The pancake disks probably fragment into fil- 
aments in which stars can subsequently form. In this case, 
HD cooling is expected to become efficient in low-density 
filaments, and then the high mass peak of the IMF would 
go down to ~ 10 Mq. 

It is found that at redshifts of 2; ~ 2, the UV back- 
ground radiation decreases with time (Irwin, McMahon, 
& Hazard 1991; Maloney 1993). The time-decreasing UV 



background radiation is likely to influence star forma- 
tion in galaxies, especially low surface brightness galax- 
ies (e.g., Ellis 1997). CorbeUi et al. (1998) studied the 
effects of the declined UV background radiation on the 
thermal evolution of the protogalaxies with low surface 
densities. They found that there is a critical redshift of 
2 ~ 1 — 2, above which the gas disks with surface densities 
10^° cm~^ < Nm < 10^^ cm~^ are gravitationally stable 
at T ~ lO''' K. Below this redshift, the declined UV radia- 
tion is shielded by the gas disks where the H2 abundance 
reaches 10~^ owing to high ionization degree by the UV 
radiation. Also, in such galaxies, the high mass peak of the 
IMF would decrease to ^ IOMq owing to the HD cooling. 

The high mass end of the IMF can influence abundance 
patterns because metal production by extremely metal- 
deficient stars is very different between lO^M© and IOMq 
(Abia et al. 2001; Umeda & Nomoto 2001; Heger & 
Woosley 2001; Schaerer 2001). For instance, the abun- 
dance pattern in the metal-poor (^ 0.05Zq) starburst 
galaxy M82 cannot be accounted for unless stars with 
> 25Mq contribute significantly to the metal enrichment 
of the galaxy (Tsuru et al. 1997; Nakamura et al. 2001; 
Umeda & Nomoto 2001). This mass scale seems to be 
consistent with the high mass peak of the present IMF 
regulated by the HD cooling. It is also consistent with the 
estimate by Hernandez & Ferrara (2001). 

As the star formation progresses, the interstellar metal- 
licity will monotonously increase with time. When the 
metallicity reaches 10"^ — 10~'^Zq, the metal cooling be- 
comes important and the thermal properties of the gas are 
changed. Thereafter, the process of star formation would 
become similar to the present-day case. In other words, 
the IMF would settle into Salpeter-hke IMF. 

We are grateful to A. Ferrara, T. Nakamoto, R. Nishi, 
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workstations at the Center for Computational Physics, 
University of Tsukuba. This work was financially sup- 
ported in part by the Grant-in- Aid for Scientific Research 
on Priority Areas of the Ministry of Education, Science, 
Sports and Culture 10147205 and 11134203 (FN). 
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Table 1 

REACTION RATE COEFFICIENTS 



Reactions 



Rate Coefficients (cm^s ^) 



Reference 



(Dl) 
(D2) 
(D3) 
(D4) 
(D5) 
(D6) 
(D7) 
(D8) 
(D9) 
(DIO) 



D+ + e ^ 
D + H+ - 
D+ + H- 
D + H^ 
D + H2 - 
HD+ + H 
D++H2 
HD + H- 
HD + H+ 
D + H+ - 



■D + hv 
^ D+ +H 
^ D + H+ 
ED + hv 
> H + HD 
^ H+ + HD 
H+ + HD 
■^Ha + D 
^ H2 + D+ 
-> HD+ + hiy 



(D12) HD+ + e ^ H + D 

(D13) D + e^B' + hp 

(D14) D++D-^2D 

(D15) H++D-^D + H 

(D16) H-+D^H + D- 

(D17) D-+H^D + H- 

(D18) D- + H ^ HD + e 



koi -- 
kD2 -- 
kos -- 
koi -- 
kD5 
koe = 
km -- 
kos - 

kD9 = 

koio 



= 3.6 X 10-i2(r/300) 
= 3.7 X 10-iorO-28exp(-43/T) 
= 3.7 X io-iotO-28 
= 1.0 X 10-25 

= 9.0 X 10-"exp(-3876/r) 
= 6.4 X 10-1" 
= 2.1 X 10-9 

= 3.2 X 10-11 exp (-3624/T) 
= 1.0 X 10-9 exp (-464/r) 
= dex[-19.38 - 1.523 logT+ 
1.118(logT)2 - 0.1269(logT)3] 

= dcx[-19.38- 1.523 logr+ 
1.118(logr)2 - 0.1269(logT)3] 

= 7.2 X io-8r-i/2 

= 3.0 X 10-i^(T/300)"-95exp(-T/9320) 

= 5.7 X 10-»(r/300)-°-5 

= 4.6 X 10-^(T/300)-°-5 

= 6.4 X 10-9(T/300)° ''i 

= 6.4 X 10-9(T/300)°-4i 

= 1.5 X 10-9(T/300)-° i 



0.75 



(Dll) D+ + H ^ HD+ + hu kon 



kD12 

kou 
koiz 
kow 
kon 
kois 



SLD98 

GP98 

GP98 

GP98 

GP98 

GP98 

GP98 

GP98 

GP98 

GP98 

GP98 

GP98 

SLD98 

SLD98 

SLD98 

SLD98 

SLD98 

SLD98 



Note. — GP98: Galli & Palla (1998); SLD98: Standi, Lepp, & Dalgarno (1998) 
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log (ric, o) 

Fig. 1. — Compaj-ison between the HD cooling rate and H2 cooling rate. The figure shows the HD-to-H2 cooling ratio in the initial gas 
density (nc,o) — initial H2 abundance (xhj.o) diagram. The HD abundance is assumed to be 1 x 10~* times the H2 abundance. The solid 
lines show the contour curves of the HD-to-H2 cooling ratio, where the equilibrium gas temperatures are determined by the condition of 
tfrag = icool aJ^6 indicated by dashed lines. If the initial parameters of the filaments are in the gray region, the HD cooling plays an 
important role in the thermal evolution of the gas (see the text for more details). 




10^ 104 10^ 



«c (cm'^) 

Fig. 2. — Evolution of (a) the temperature, (b) H2 {solid line) and HD {dashed line) abundances, (c) total {solid line), H2 {dashed line), and 
HD {dotted line) cooling rates, (d) the cooling time {solid line), the contraction time {dashed line), and the fragmentation time {dotted line), 
respectively, as a function of the central density [ric = Pc/{tJ''mH}] for the model with (nc,o, To, /, xhj^o) = (lOcm^*^, 200 K, 2, 3 x 10^'^). 
Note that the definition of the fragmentation time is 2.5 times larger than that in Paper H. 
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/ie,o(cmO /ie,o(cmO 

Fig. 5. — The distributions of fragment mass for the models with (a) Tq = 300 K and xhj.o = 1 x 10-3 and (b) Tq = 200 K and 
xh2,o = 3 X IQ-^. (To is close to the equilibrium temperature for each value of xu^fi-) The abscissa and ordinate denote the initial central 
density and the parameter /, respectively. The solid lines denote the contours of the fragment mass which are labeled with adjacent numbers 
in units of Mq. A thick dashed line shows the lines at which Ahd = the epoch of fragmentation. In the left regions of the dashed 

lines, HD cooling is more efficient than H2 cooling. The dot-dashed lines show the line masses of the filaments formed by fragmentation of 
self-gravitating gas disks (see the text for more details), where /ieq is a line mass of an equilibrium filament. The number attax;hed by each 
line denotes the line mass in units of Heq- 



